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Abstract 

Coulomb wave functions are difficult to compute numerically for extremely low energies, even with direct 
numerical integration. Hence, it is more convenient to use asymptotic formulas in this region. It is the object 
of this paper to derive analytical asymptotic formulas valid for arbitrary energies and partial waves. Moreover, 
it is possible to extend these formulas for complex values of parameters. 

PACS: 02.30.Fn, 02.30.Gp, 03.65. Ge 

1 Introduction 

Coulomb wave functions are the typical example of functions being both analytic and very difficult to compute. 
They are part of the rare cases for which second-order differential equations can be analytically solved and are 
expressed with confluent hypergeometric functions [1]: 
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where I is the angular momentum of the wave function, r\ its Sommerfeld parameter, u> can be equal to ±1 
in eas. (|1I2|) . and where the constants of normalization Ci(rj) (Gamow factor) and crt{r]) (Coulomb phase shift) 
appear pQ. All computational difficulty arises from the presence of confluent hypergeometric functions \F\ and 
2-Fo- They indeed vary by many orders of magnitude for smooth variations of parameters and are moreover 
subject to a cut in the complex plane when analytically continued [T]. Codes handling arbitrary complex 
parameters in Coulomb wave function computation have been published in both Fortran [5] and more recently in 
C++ |3J languages. While the former uses only analytical methods such as power series and continued fractions, 
the latter included direct integration as well, which has considerably extended the numerical domain of definition 
for which implementation of Coulomb wave functions is stable [3]. However, even with this amelioration, results 
become unreliable for very large values of |Im(^)| and/or \r/\ [5J. 

In ref.[l], a uniform approximation for Coulomb wave functions has been presented, which, however, demands 
I = 0. This prevents partial decomposition of wave functions, useful for reaction cross section calculation 
[5]. Moreover, approximations for which I is arbitrary would be of interest for relativistic calculations, where 
Coulomb wave functions are expressed also with the confluent hypergeometric functions appearing in eas. (|H2p 
®- 



2 Uniform approximation of Coulomb wave functions 

In order to alleviate instabilities encountered in Coulomb wave function implementation, we will derive analytic 
formulas valid for large ?y, £ > and p > firstly, which will be secondly analytically continued to complex 



1 



values. The standard method therein is to use uniform WKB approximation, where divergences occurring at 
the turning point, denoted in the following as pt, are removed through the use of Airy functions [7]. While the 
method described in ref.JTj demands in the general case to deal with non-analytical integrals, those appearing 
for the Coulomb problem can be calculated exactly with elementary functions. To apply this method, one firstly 
writes the exact ansatz verified by Coulomb wave functions, from which approximations can be effected: 
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where Ai(z) and Bi(z) are the standard regular and irregular Airy functions [T] (their normalization will be 
justified afterward, as well as the reason why Fg v (Ge n ) possesses no Bi (Ai) component), and where <j>{x) verifies 
the following third-order non-linear differential equation: 
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where the parameter a has been introduced for convenience. For large tj, p t — > +oo, so that terms proportional 
to p^ 2 in eq. (fT0|) can be neglected. Hence, tf>(x) verifies asymptotically a non-linear first order equation. 
Elimination of turning point divergence moreover demands that 0(0) = 0, so that its approximate equation and 
solution read for real x: 

^4>{x)i = I J r f T + u T,^, dt,x>0, 
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This approximation has been considered in ref.[J] in the particular case I = 0, implying a = 0, for which 
4>{x) is equal to a particularly simple expression. For x — ► +oo, one can verify easily that (2/3)0(a;) 3 / 2 ~ x, 
which justifies the ansatz of cqs. (|6I7|1 . where they reduce to standard sine-cosine approximation pQ. If a ^ 0, it 
happens that eq. (|13[) not only can be integrated analytically, but also in a very concise way: 
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One can check that the I = case described in ref.[l] is properly obtained using a = in eq. (fl"4"|) . 4>'{x) is 
obtained by eq. p^]) and the condition that <p'(x) > 0, immediate from eq. p3|) . Note that for x ~ 0, eq. (fT4)l 
becomes numerically unstable, so that it is preferrable to use the asymptotic formulas verified by cf>(x) and 
4>'(x) therein, which are (f>(x) ~ (1 + a) 1 / 3 ^ and 4>'{x) ~ (1 + a) 1 / 3 . FL{p) and G' e (p) are obtained by a simple 
differentiation of cqs. (|6!7p . It has been noticed numerically, however, that their term proportional to <f>"(x) 
is not negligible and should be kept in approximate formula. (f)"(x) is obtained by a simple differentiation of 
eq.(P]). 
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Extension of eq. (fT4")) to complex arguments demands caution. On the one hand, it is usually sufficient therein 
to replace the conditions x > and x < by Re(ir) > and Re(x) < respectively, especially if considered 
imaginary parts are small in modulus (see also ref. [5] for computational methods of the Airy function in the 
complex plane). On the other hand, however, one has to pay attention to the different cuts obeyed by both 
Coulomb wave functions and the elementary functions of eq. lfTI")) . The theoretical behavior of Coulomb wave 
functions in the vicinity of their cut has been studied in ref. [9] . The simplest method to avoid problems generated 
by cuts is to consider complex contours which never cross the negative real axis, so that Coulomb wave functions 
are continuous therein. Then, one just has to modify formulas of eq.([rjj so that they are continuous on these 
contours if cuts of elementary functions therein are encountered. They are straightforward to treat, as cuts 
appear only by way of log, inverse circular/hyperbolic and power functions. 

3 Numerical examples 

In order to show the efficiency of the approximation presented in eq. (fTi|) , one will consider both sets of param- 
eters. The first set consists in the real values 1 = 1 and r\ = 10, for which p > 0. The second set reads I = 2 + i 
and 77 = 10 + i, with p = \p\e l7T ^ 4 . No cut problem appears for these complex values in eq. (fT4")) . |Re[F^(p)]|, 
|Re[G&j(p)]| and analog values related to Coulomb wave functions derivatives are depicted in log scale, in fig.(H]) 
for the real set of parameters and in fig.© for the complex set of parameters. They are compared to the exact 
functions calculated numerically with the code of ref. [3] , with which Airy functions are calculated as well as 
linear combinations of Bessel functions pQ. It is clear from these figures that the approximation provided by 
eq. (|14p is very good, even though i] is not very large. The relative error of the approximate formula of eq. (|14p 
is ~1% for most p values, which was expected as precision of the approximation is of the order of \pt\ 2 (see 
eq.CHj). 

4 Conclusion 

The uniform Coulomb wave function approximation of ref. 4 , valid for r\ — > +00 and t = 0, has been generalized 
to both arbitrary partial wave and complex parameters, and has been checked numerically to be reliable. 
It provides a useful alternative to exact computation of low-energy Coulomb wave functions, which can be 
numerically costly or unstable. In particular, partial wave decomposition can be effected with the proposed 
uniform approximation of low-energy wave functions. 
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Figure 1: Absolute values of Fi v (p), Gi v (p), F' lr) (p) and G'^Jp) for I — 2 and 77 = 10. Exact calculation provided 
by the code of ref.[3] and approximation issued from eq. (|14|) are indistinguishable. 
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Figure 2: Absolute values of Re[F er ,(p)], Re[Gi v (p)}, Re[F^(p)] and Re[G' ir] (p)} for £ = 2 + i, r) = 10 + i and 
arg(p) = 7r/4. Exact calculation provided by the code of ref.[3] is provided as straight lines and approximation 
issued from eq. (fTl| as dashed lines. Differences are visible only for |Re[G^(p)]| and |Re[GL(p)]|. 
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